Macroscopicity of Mechanical Quantum Superposition States 
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We propose an experimentally accessible, objective measure for the macroscopicity of superposi- 
tion states in mechanical quantum systems. Based on the observable consequences of a minimal, 
macrorealist extension of quantum mechanics, it allows one to quantify the degree of macroscopicity 
achieved in different experiments. 



Introduction. — Experiments probing the quantum su- 
perposition principle at the borderline to classical me- 
chanics are a driving force of modern physics. This 
includes the demonstration of superposition states of 
counterrunning currents involving 10 14 electrons 0, 
of Bose-Einstein condensed atoms Q, and complex 
molecules 

Various measures have been suggested for the size of 
superposition states involving macroscopically distinct 
properties of complex quantum systems [a-Hlj. Most of 
them refer to specific types or representations of quantum 
states, or count the operational resources required to an- 
alyze them. While most proposals seem to be grounded 
in a common information-theoretic framework [l2j |. we 
still lack a method of attributing a definite and unbiased 
measure to all experimental tests of the quantum super- 
position principle. 

The task to define a macroscopicity measure 'within' 
quantum theory is confounded by a fundamental prob- 
lem: We are free to decompose a many particle Hilbert 
space into different tensor products, such that a com- 
plicated single-particle representation of a wave function 
may look mundane after a change of variables to collec- 
tive degrees of freedom. This highlights the problems of 
an ad hoc selection of distinguished observables. 

In view of this, we propose to call a quantum state 
of a mechanical system the more macroscopic the bet- 
ter its experimental demonstration allows one to rule 
out even a minimal modification of quantum mechanics, 
which would predict a failure of the superposition princi- 
ple on the macroscale. Turning this characterization into 
a definite measure requires one to specify the minimal 
modification. 

Fortunately, it is not necessary to worry about the de- 
tails of possible nonlinear or stochastic additions to the 
Schrodinger equation, which might embody the coarse- 
grained effects of a deeper theory, say, incorporating 
gravitation or a granular space-time I13-ll5l . or might 
represent a fundamental stochasticity 16-18]. All that 
matters empirically are their observable consequences, 
described by the dynamics of the many-body density 
operator. We argue that basic consistency, symmetry, 
and scaling arguments lead to an explicit, parametriz- 
able characterization of the impact of a minimally in- 
trusive modification. The fact that no evidence of such 



physics 'beyond the Schrodinger equation' is seen in a 
quantum experiment rules out a certain parameter re- 
gion. For a superposition state in a different experiment 
to be more macroscopic, its demonstration must exclude 
a larger parameter region, implying that possible modi- 
fications must be even weaker. Diverse experiments can 
thus be compared without prejudice. 

Minimal modification of quantum mechanics. — The 
modification must serve to 'classicalize' the state evolu- 
tion in the sense that superpositions of macroscopically 
distinct mechanical states are turned rapidly into mix- 
tures. The operational description of quantum theory, 
based on the state operator p, its completely positive and 
trace-preserving time evolution, and a consistent rule of 
assigning probabilities to measurements (l9) . allows one 
to treat (nonrelativistic) quantum and classical mechan- 
ics in a common general formalism. It is therefore natural 
to account for an objective modification of the quantum 
time evolution in the framework of dynamical semigroups 
20]. That is, the effect of the modification can be ex- 
pressed as a generator Cn added to the von Neumann 
equation for the state of motion pj^ of an arbitrary sys- 
tem of N particles, dtPN — [H, pj^]/ih 4- CnPn- 

In addition, we require the modification (i) to be in- 
variant under Galilean transformations, avoiding a dis- 
tinguished frame of reference, (ii) to leave the exchange 
symmetry of identical particles unaffected, (iii) to respect 
the 'innocent bystander' condition that adding an uncor- 
related system leaves the reduced state unchanged, and 
(iv) to display scale invariance with respect to the center- 
of-mass (cm.) of a compound system. We will see that 
these requirements essentially determine the form of a 
possible minimal modification. 

Let us first consider an elementary particle of mass m 
using the formalism of quantum dynamical semigroups. 
A theorem by Holevo [21| states that any Galilean invari- 
ant addition to the von Neumann equation for the state 
of motion p must have the form dtp = [H, p]/ih + C\p 
with 
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if one disregards unbounded diffusion terms, which would 
yield a substantially more drastic modification. The op- 
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effect a translation s and a velocity boost v of the ele- 
mentary particle, while g (s, q) is a positive, isotropic and 
normalized phase-space distribution, whose standard de- 
viations for the position and the momentum variable will 
be denoted by o~ s and o~ q . The von Neumann equation is 
reobtaincd for a s = a q = 0. 

The modification (jTJ serves its purpose of classicalizing 
the motion of a single particle: It effects a decay of the 
position and the momentum off-diagonal matrix elements 
of p. The parameter r provides the corresponding time 
scale for those matrix elements which are more than the 
critical length scale H/a q off the diagonal in position, or 
more than h/o~ s off in momentum, whereas smaller-scale 
coherences may survive much longer. Delocalized super- 
position states thus get localized in phase space as time 
evolves, ultimately rendering the phase-space represen- 
tation of p indistinguishable from an equivalent classical 
distribution. At the same time, the modification (TTJ) in- 
duces a position and momentum diffusion, implying that 
any bound particle gradually gains energy. For harmonic 
binding potentials with frequency uj the energy increases 
as a q /2m + muj 2 a1l2 per unit of time r. 

Many-particle description. — The requirement of 
Galilean invariance in a general mechanical system of N 
particles implies that the phase-space translation opera- 
tors must effect a net shift of the center-of-mass coordi- 
nates by s and v. On the other hand, the scale invariance 
conditions with respect to an innocent bystander (iii) and 
to the center-of-mass (iv) require that the equation for 
the iV-particle state reduces to the single-particle form 
whether one traces over the other N — 1 particles or over 
the relative coordinates in a compound object of rigidly 
bound contsituents; in the latter case the single mass 
should be replaced by the total mass M = J2 n m n- This 
is achieved by composing the TV-particle operators as the 
weighted sum of the single-particle operators ©, 
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where m c is an arbitrary reference mass; see the Ap- 
pendix for details. 

We note that the operators ([3]) conserve the exchange 
symmetry of a quantum state. The corresponding N- 
particle equation 
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thus leaves boson and fermion statistics invariant (ii). 
The equation is completely determined once we specify 
the mass m e , the coherence time parameter r e , and the 



normalized distribution function g c (s, q) for the reference 
particle. The innocent bystander condition (iii) guaran- 
tees that no correlations are introduced between different 
(possibly uncorrelated or even far apart) subsets of parti- 
cles. Property (iv), on the other hand, admits the single- 
particle description (TTJ not only for elementary point par- 
ticles (e.g. electrons) but also for compound objects such 
as atoms, molecules or even solids. 

The classicalization of the center-of-mass motion of an 
extended compound object of total mass M can be ap- 
proximated by the single-particle form (TTJ) , if the relative 
motion of the constituents around their rigidly bound 
equilibrium positions can be neglected. The Fourier 
transform g (q) — J d 3 x g (x) e~ tq ' x / h of the mass den- 
sity g (x) of the compound modifies the rate and the 
phase-space distribution of the effective center-of-mass 
classicalization, 
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The effective coherence time r depends on the relation 
between the size of the compound and the critical length 
scale H/a q of the reference distribution g c ; the effective 
distribution g remains normalized. The description fails 
as soon as the relative motion of the constituents must 
be taken into account. The inner structure of nuclei 
becomes relevant for femtometer-scale distribution func- 
tions, h/o~ q < lOfm < (rn/m e )a s , with m « 1 amu 
the nucleon mass. Here ends the domain of nonrela- 
tivistic quantum mechanics, and with it the validity of 
our approach. We thus restrict its parameters to about 
o~ s < 20 pm and h/a q > 10 fm, noting that the macro- 
scopicities will not change if these bounds are varied by 
a few orders of magnitude. 

The restriction to a single reference distribution 
g e (s,q) in (fj|, as opposed to individual distributions 
for different types of particles, yields a universal single- 
particle description (fTJ). The choice of the reference 
mass m c is arbitrary, since the coherence time param- 
eter and the distribution rescale to r = t c (m c /m) 2 and 
g(s,q) = (m/m e ) g e (ms/m e ,q) for a point particle of 
different mass m, as follows from This renders the 
translation s negligible for heavy objects, m 3> m e . 

In the following we will use the electron as reference 
particle fixing both r c and g c . Moreover, we take g c to 
be a Gaussian distribution in s and q, fully specified by 
the standard deviations <r s and <r q . The latter determine 
the main behavior of the classicalization effect; a more 
involved description with additional parameters would 
complicate matters without significantly modifying the 
generic behavior. 

It is remarkable that a special form of Eq. Q , which 
we arrived at using the assumptions (i)-(iv), describes 
the observable consequences of the theory of "continu- 
ous spontaneous localization" (CSL) [22[ if one takes 
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FIG. 1. Lower bounds on the time parameter r c , as set by 
various experiments. The calculations are done for the rele- 
vant range of critical length scales h/a q , and at a s — 20pm. 
The solid line corresponds to the atom interferometer of [2q |: 
it rules out all time parameters r e below the curve. Future 
experiments may exclude a larger set e.g. by interference of 
10 5 -10 7 amu gold clusters [2^] (dashed lines) or of micromirror 
motion [301 ] (dash-dotted line). The dotted line corresponds 
to demonstrated persistent current superpostions in a SQUID 
loop [|. The shaded region represents the excluded r e by a 
conceivable classical measurement of less than 1 /xK/s tem- 
perature increase in a Rb gas. 



cr s = [23(. This shows that one can set up explicit 
theories which modify the dynamics on the level of the 
Schrodinger equation and whose observable consequences 
fit into the present framework |24l-[27|. The stochastic 
Schrodinger equation in [l|| |22|| may thus be seen as 
one example, but not the most general form, of a the- 
ory which yields a minimal modification in the sense de- 
scribed above. 

Assessing superposition states. — The experimental 
demonstration of quantum coherence in a mechanical de- 
gree of freedom rules out a certain parameter region of 
the classicalizing modification, i.e. it provides a lower 
bound of the time parameter r c for any fixed value of 
a s and o~ q . For a superposition state in a different exper- 
iment to be more macroscopic, its demonstration must 
exclude a larger set of r e , implying that the modification 
must be even weaker. 

Figure Q] shows the greatest excluded r e for a number 
of different setups, as a function of the critical length 
scale hjdq and at fixed a s = 20 pm. The solid and the 
dotted curve correspond to exemplary modern matter- 
wave experiments: The interference of cesium atoms in 
free fall over hundreds of milliseconds (solid line) [28j . 
and the superposition of counterpropagating currents of 
10 14 superconducting electrons in a Josephson ring (dot- 
ted line) [l|. The dashed and dash-dotted lines illus- 
trate what would be achieved in proposed superposition 
experiments with nanoclusters [23 | or micromirrors [30| . 
Whereas the value of a s matters for the SQUID experi- 



ment, it is not important for the other cases due to the 
large masses involved. Our results then resemble the pre- 
dictions of a CSL model with varying localization length 
h/o~ q . Detailed results on each experiment are reported 
in the Appendix. 

One observes a common feature of all quantum curves 
inFig.[TJ they saturate or assume a local maximum. This 
is because the classicalizing modification (j4]) is bounded 
in the operator norm, and any given position or momen- 
tum superposition state of a total mass M thus survives 
at least for a time r (m c /M) 2 . 

The interferometer results (solid and dashed lines) 
reach the maximum where the critical length scale H/o~ q 
is comparable to the interference path separation. This 
is where the solid line saturates, in accordance with the 
CSL results in (2(| . The dashed line drops at length scales 
smaller than the size of the interfering object, when only 
a fraction of its mass contributes to its center-of-mass co- 
herence time, as given by Eq. ([5]). The latter also holds 
if the object is larger than the path separation (dash- 
dotted curve). For smaller values of h/a q the slope is 
mainly determined by the mass density g(x) of the in- 
terfering object, while in the diffusive limit of large values 
it solely depends on <j q . 

The superposition of persistent currents probed in the 
SQUID experiment [ij can be described by displaced 
Fermi spheres of Cooper-paired electrons |3l| . The clas- 
sicalization gradually redistributes and dephases elec- 
trons between the Fermi spheres, thereby undermining 
the quantum coherence (see Appendix). At large mo- 
mentum spreads a q the effect is governed by the redistri- 
bution of electrons, as would be the case in spontaneous 
localization [32j . The dotted curve assumes its maximum 
at a value o~ q where the redistribution covers all electrons 
in the Fermi sphere. This effect vanishes for smaller a q 
once the superconducting energy gap can no longer be 
overcome; in this limit the classicalization effect is gov- 
erned by the dephasing that is induced by the position 
diffusion with spread o~ s . 

Because of the diffusion effect inherent in the clas- 
sicalizing modification there are also 'classical' experi- 
ments not explicitly demonstrating quantum behavior 
that can be used to narrow down the range of plau- 
sible coherence time parameters. This is indicated by 
the shaded area in Fig. Q] which represents the values of 
t~q excluded by an anticipated precision measurement of 
the temperature increase of a dilute gas of Rb atoms, 
0- 2 /2mRbT Rb < 1.5fcs x 1/iK/s. 

Macroscopicity measure. — In view of the parameter 
bounds displayed by Fig. Q] we suggest to quantify the 
macroscopicity of a superposition state realized in an 
experiment by the greatest excluded time parameter r e 
of the modification respecting the above parameter 
bounds. To set a scale, we take the logarithm of r in 
units of seconds as the measure of macroscopicity, 
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FIG. 2. Timeline of macroscopicities reached in quantum 
superposition experiments. (See the Appendix for details.) 
The squares, the triangles and the dots represent interfer- 



ence experiments with neutrons 33, 34], atoms [2 
atom BECs 3], and molecules \4, 139I444T ]. respectively. One 
notes that BECs do not substantially exceed the macroscop- 
icities achieved with atom interferometers. This is due to the 
single-particle nature of the condensate wave function. The 
many-particle state is more involved in the case of superpo- 
sition experiments with persistent supercurrent states in a 
large SQUID loop [3, @| , as represented by the stars. How- 
ever, despite the large number of Cooper pairs contributing to 
the current superpositions in SQUIDs, such experiments lag 
behind in macroscopicity due to the small coherence times 
observed. 



That is, a positive value for /i is obtained if one demon- 
strates an electron to behave like a wave for more than 
one second, or of a proton for about a microsecond. 

A simple approximate expression for the macroscopic- 
ity /i is obtained for interference experiments with point 
particles or with compound objects of total mass M, 
whose size is much smaller than the path separation. The 
single-particle modification (JXJ) predicts an exponential 

decay of coherence with time scale r = r e (m e /M) , a 
mass dependence also obtained in the CSL case |24j - t26| . 
This is to be compared with the period t during which co- 
herence is maintained in the experiment. Measuring with 
confidence a fraction / < 1 of the expected interference 
visibility, one gets 
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That is to say, if one measures 30% contrast in an in- 
terference experiment where the visibility is predicted to 
be, say, 60%, then / = 0.5 must be used in the above 
expression. 

Macroscopicity of specific experiments. — In Fig. [2] 
we present the macroscopicities attained in a selec- 
tion of quantum experiments versus their publication 
date. They include tests of the superposition principle 
with neutrons, electrons, individual and Bose-condensed 
atoms, and molecules. Details on the calculations for 
specific experiments can be found in the Appendix. 



Oscillating micromembrane 
Hypothetical large SQUID 
Talbot-Lau interference [291 at 10 5 amu 
Satellite atom (Cs) interferometer [45] 
Oscillating micromirror [30J 
Nanosphere interference [46] 
Talbot-Lau interference [29] at 10 8 amu 
Schrodinger gedanken experiment 



11.5 
14.5 
14.5 
14.5 
19.0 
20.5 
23.3 
~ 57 



TABLE I. Expected macroscopicities for various proposed and 
hypothetical quantum superposition experiments. See the Ap- 
pendix for details. The oscillating micromembrane setup [43] 
will reach the stated /Li-value if coherence between the zero- 
and one-phonon state can be observed for over 1000 oscilla- 
tion cycles. For the SQUID experiment we assume a loop 
length of 20 mm, a wire cross section of 100 pm 2 , and 1ms 
coherence time. In the gedanken experiment an idealized cat 
of 4 kg is kept in a spatial superposition of 10 cm distance for 
Is. 



State-of-the-art interferometers achieve macroscopic- 
ities of up to ji « 12, and various ideas to surpass 
this value with future experiments have been suggested. 
As can be seen from Tab. HI the most promising pro- 
posals from the perspective of the macroscopicity mea- 
sure employ oscillating micromirrors 30] and nanoclus- 
ters [29. 46]. Their huge mass would trump a conceivable 
SQUID experiment with more than 10 17 electrons or an 
atom interferometer hovering in free space with an inter- 
rogation time of one hour [45[ . 

Nevertheless, there are more than 30 orders of magni- 
tude between experiments conceivable with present-day 
technology (/i = 12-24) and something as manifestly 
macroscopic as an ordinary house cat (p ~ 57). 

Conclusion. — Using the measure proposed in this ar- 
ticle any experiment testing the superposition princi- 
ple in mechanical degrees of freedom can be quantified 
and compared. By definition it answers an empirically 
relevant question, namely to what extent an observa- 
tion serves to exclude minimally invasive modifications 
of quantum mechanics that produce classical behavior 
on the macroscale. As such, the measure follows di- 
rectly from basic symmetry and consistency arguments, 
and confers physical meaning on the abstract notion of 
macroscopicity of a quantum system. 

The proposed measure does not depend on how a com- 
pound mechanical object is decomposed into elementary 
mass units. For instance, an interfering fullerene bucky- 
ball might be described in terms of 60 carbon atoms or 
equally of 1080 nucleons and electrons, and both descrip- 
tions should consistently lead to the same macroscopicity 
value for the overall state of the molecule. This issue, 
which was not addressed in previous studies, is explic- 
itly taken into account in our approach. Moreover, we 
do not refer to specific classes of quantum states, or to 
preferred measurement operations or observables, render- 
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ing the measure of macroscopicity applicable to arbitrary 
mechanical systems. 

The last 20 years have witnessed a remarkable rise in 
demonstrated macroscopicities. Yet, new experimental 
strategies for quantum tests, in particular using nan- 
oclusters and microresonators, may soon venture deeper 
into the macroworld. As more and more effort is put into 



this field, we may well experience an unprecedented leap 
towards the macroscopic domain. 
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APPENDIX 



Here we provide a more extensive motivation for the many-body form (5) of the modification, and a justification 
for our estimates of the macroscopicities reached in the various experiments. We note that a more detailed analysis 
of the individual experimental setups might yield slightly different values, but would not change the overall picture. 



Appendix A: Many-particle form of the classicalizing modification 



The requirement of invariance under Galilean transformations (i) implies that the generator describing the effect 
of a minimal modification is of the form 



C N p = \d 3 s& 3 v 
The operators L(s,v) must satisfy 



L(a,v)p\J (s,v) - - {0 (s,v)L(s,v),p} 



exp 



(P • s' - Mv' ■ X) 



L (s, v) exp 



(P ■ s' - Mv' ■ X) 



exp 



h 



(v ■ s' 



l(s,v) 



(Al) 



(A2) 



with M = rn n the total mass, to an arbitrary reference mass, and X, P the center-of-mass position and 

momentum operators. By switching to center-of-mass and relative coordinates it follows from Eq. (|A2j) that the 
L(s,v) induce a net shift of the center-of-mass position and momentum by m c s/M and m e v 7 respectively. However, 
it remains unspecified how the net shift is to be distributed amongst the N constituents of the system. This freedom 
is constrained by the additional assumptions (ii)-(iv). 

Assumption (iii) means that the A-particle form (|A1[) must always reduce to the single-particle form for the nth 
particle, Eq. (1) in the main text, if one traces over the other N — 1 constituents, trjv-i (Cnp) = £itrjv-i (p). That is 
to say, we assign to each particle species of mass m n an individual time parameter r n and a positive, normalized and 
isotropic distribution function g n (s,q). They constitute the free parameters of the single-particle form, as discussed 
in the main text. At the other end, assumption (iv) recovers the single-particle form (1) for the center-of-mass degree 
of freedom provided that the N constituents are well localized at fixed equilibrium positions close to the center. We 
denote the corresponding classicalization parameters by and g( N \ Finally, the TV-particle operators L (s, v) must 
be symmetrized in the case of indistinguishable particles (ii) . 

As a first guess one might think that the L (s, v) should be proportional to the unitary A-particle Weyl operators, 
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(A3) 



This way the phase-space shift would be distributed equally among all participating particles. The assumptions (ii) 
and (iv) would be fulfilled by construction and, due to the prefactor in (|A3|) . one would reobtain the center-of-mass 
distribution g^ N > (s,q) in (IA1|) . However, the operators (|A3[) would leave the relative motion of any constituent 
subsystem entirely unaffected, irrespectively of the overall size and extension of the iV-body system. Moreover, 
assumption (iii) is met only if = r„ and (s,q) — (m n /M) 3 g n (s,m n q/M) for all n, so that the effective 
classicalization rate I/t'"' would not increase with the system size. Therefore, the operators (|A3[) cannot induce 
classical behavior at the macroscale leaving at the same time microscopic systems unaffected. 
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Rather than dividing the phase-space shift among many particles, one may as well compose a solution of 
single-particle translations, 
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(A4) 



The sign (±) may differ for distinguishable particles, and it may also depend on s and v. These operators fulfill (ii) 
and (iii) by definition. Moreover, in the case of a compact compound, where x„ rs X and p n ~ m n P/M, we recover 
condition (iv), with the center-of-mass parameters determined by 



N 



-9n 



n=l 



rn. 



s,q 



(A5) 



Each constituent contributes to the collective classicalization of the center-of-mass variables. For instance, in the 
case of N indistinguishable particles we find that the classicalization rate is amplified by = N 2 /t. A formal 

proof of the operator expression (IA4|) can be obtained in the picture of second quantization if one takes the operators 
L (s, v) to be a combination of single-particle terms, i.e. a bilinear form in the annihilation and the creation operators 
of a given particle species. The ./V 2 -scaling of the classicalization rate then follows immediately. 

In general the sign factor (±) n and the single-particle distributions g n might still differ from species to species. This 
is where the scale-invariance argument (iv) can be invoked once again: If we allowed for different signs and distribution 
functions for different point-particle species, we would have to define a fixed set of reference point particles, i.e. single 
out a distinguished many-body representation of any composite mechanical system. Moreover, we would end up with 
different descriptions of the modified time evolution of approximate point-like compound particles depending on their 
composition. 

In order to guarantee true scale invariance, by avoiding an ambiguous treatment of approximate point particles, 
one must therefore relate the extensive nature of the modification in (|A5|) to the elementary extensive property of 
mechanical systems: their mass. We notice that the summands in (| A5|) contribute with a priori different mass scales 
m n , each of which might be composed of further sub-units of mass, and so on. A unified description for any mass 
scale is obtained only by introducing a single reference time parameter r e and distribution function g c associated with 
a fixed reference mass m c . The composition rule (|A5[) then holds naturally once we identify 



(±), 



I mi ( m c \ m n / 1 . . 
—3 — 9n s,q = \ —9e{s,q). 



(A6) 



The classicalization effect now scales uniformly with mass, irrespectively of the types of particles involved. 



Appendix B: Diffraction at gratings and double-slits 



The macroscopicity observed in the matter- wave diffraction experiments at gratings and double-slits [33H35L 1371 . 1411 . 
HH can be estimated from the height of the first order diffraction peak in the recorded signals. Given the transmission 
function t(x) of a one-dimensional iV-slit-grating with slit distance d and opening width iv, we find the interference 
signal of a monochromatic point source (velocity v z ) in the paraxial approximation 



I Vz (x) oc / dxidx2 R (xi — xi) t{x\)t* {x2) exp 



2h 



x\ 1x{x\ — x%) 



T 



T-2 



(Bl) 



The time T = T\Tij^T\ + T2) is determined by the times-of-flight T\ from the source to the grating and T2 from the 
grating to the screen. In a horizontal alignment they are related to the respective distances Li 2 = w z T 12 . The signal 
(|BI I must be averaged with respect to the distributions of the velocities v z , and over the extensions S and D of the 
source slit and the detector. The single-particle classicalization as described by equation (I) in the main text blurs 
the interference signal by the factor 



R(x) = exp 



Ti 



9 id xz 



Ti 



- 1 



T 



9iD xz, 



T-2 



1 



(B2) 



involving the reduced Fourier transform ^yd{x,p) — g(xe x ,pe x ) of the distribution function g. The unperturbed 
fringe pattern exhibits diffraction maxima at screen coordinates close to integer multiples of x = hT2/md. The clas- 
sicalization affects them strongest if the contributing interference paths are completely resolved by the critical length 
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scale K/(Tq. In this limit a flat background is added to the overall signal, and the longitudinal velocity distribution 
affects only weakly the reduction of the first diffraction maximum; the latter can thus be used to extract the macro- 
scopicity /x for each experiment according to equation (9) of the main text. The required parameter / is estimated 
by the ratio of the measured height of the diffraction peak and its unperturbed theoretical value, both normalized 
to the integrated signal. The following table contains all required parameters of the different experiments, including 
references to the data used for the comparison. 



Ref. 


Fig. 


Li/m 


1/2 /m 


N 


d/ fim 


w/d S / /im 


D/fim 


<M/~ 


Av z /(v z ) 


m/amu 


/ 




[33] 


9 


4.0 


5.7 


2 


107 




10 


30 


907 




1 


0.6 


4.8 


m 


7 


5.0 


5.0 


2 


126 


0.17 


15 


30 


216 


0.05 


1 


0.9 


6.2 


[35] 


2b 


1.0 


1.5 


50 


0.2 


0.5 


10 


25 


1000 


0.12 


23 


0.5 


6.8 


[48] 


3b 


0.08 


0.11 


2 


6 


0.33 


20 


20 






20 


0.8 


9.1 


[37J 


1 


0.45 


0.52 


100 


0.1 


0.43 


10 


25 


396 


0.1 


84 


0.8 


8.3 


[41j 


2a 


1.14 


1.25 


100 


0.1 


0.38 


10 


8 


226 


0.6 


720 


0.6 


10.6 



The neutron interference at a biprism observed in [331 ] can be related to the coherent superposition of two virtual 
sources separated by d = 107 /Ltm, 9.7 m away from the detector. The resulting fringe pattern thus resembles a double- 
slit pattern. The authors of [33[ present the measured data and a fitted theory curve in Fig. 9. The data deviates 
from the predicted height of 200 a.u. of the first diffraction order by roughly 50 a.u. Subtracting a dark count rate of 
60 a. u. leads to the estimate / ~ 1 — 50/140 ~ 0.6. The experiment [48| is a vertically aligned interferometer, where 
neon atoms are released from a trap. They fall through a double-slit and into a detector within t ss 200 ms. We 
obtain f by comparing the measured diffraction peak in Fig. 3b with the theoretical model of [Hf in Fig. 3g. For 
[3314351 T37ll48l we use a Gaussian velocity distribution with the mean (v z ) and the FWHM Av z , as specified in the 
table. For [41[ we use the distribution provided in the article, and we account for the special detection scheme by 
replacing the detector slit by a Gaussian laser focus of waist D. The dispersive interaction between the particles and 
the grating walls is taken into account for [37], EJ by reducing the effective slit opening size. 

The macroscopicity of the proposed optical double-slit experiment [46j with silica nanospheres is estimated by 
considering the Fourier amplitude which corresponds to the expected double-slit fringe oscillation (|Blj) . The classical- 
ization modifies it by the factor R (d) , with the largest proposed value for the slit distance d = 52 nm. We evaluate the 
macroscopicity by modeling the particles as homogeneous spheres (g = 2200 kg/m 3 ) of 20 nm radius (see equations 
(6) and (7) of the main text), and by assuming that at least 50% of the fringe amplitude is observed. 

Talbot-Lau interference with molecules and clusters pi. [29T I42TI44I ] can be treated in a similar manner. The sinusoidal 
fringe visibility V s i n of a symmetric setup (Ti = T% = T) is reduced to R (hT/md) V S i n . Judging from the error bars 
at high visibilities, we assume that the measurements are compatible with at least 90% of the prediction for C70 
molecules (Fi& 3 in (H), 90% for C 6 oF 4 8 (Fig. 5 in H|), and 80 % for PFNS8 (Fig. 4b in Q). The fringe pattern 
observed in [43[ with C60F48 (Fig. 4) corresponds to 75% of the predicted visibility. 



Appendix C: Ramsey-Borde interference with I2 molecules 

In the experiment (39l | a beam of I2 molecules (m = 254 amu) passes two pairs of counterpropagating running-wave 
laser beams, as described in detail in (4t|. Two paths through the setup contribute to the recorded Ramsey fringe 
pattern, as shown in Fig. 2 of Ref. (39[. We must include, however, a significant contribution from two further paths 
to the signal [49j. Their interference is washed out over the transverse velocity distribution of the molecule beam. 
We assume that they add an offset to the most pronounced central fringe of the (*)-curve in Fig. 2 of [3^|. If all 
four paths contribute by roughly the same weight we must halve the offset of the central fringe in the diagram, which 
yields a two-path fringe visibility of / « 400 a.u./(2400 — 1200) a.u. = 0.33. The passage time is determined by the 
(35 + 2) mm length of the interferometer and the mean molecular velocity of 350 m/s. This yields /1 = 7.3, according 
to equation (9) of the main text. 



Appendix D: Mach-Zehnder-type interference 

The two atom interferometers featuring the greatest macroscopicity in Fig. 2 [2^, |38| . as well as the proposed 
satellite atom interferometer [45| listed in Tab. 1, are optical Mach-Zehnder-type geometries, which could in principle 
yield close to 100% fringe contrast. We use the recorded fringe visibilities / = 0.62 (Fig. 19 in [28]) and / = 0.33 
(Fig. 3 in [38]), and a hypothetical value of / = 0.5 for the proposal [45[. In all three cases the interfering particles 
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are 133 Cs atoms, and the interrogation time is given by twice the pulse separation time T. The respective values are 
T = 160 ms, 400 ms, and 2000 s. 

In the Na2 molecule interferometer [!(| the Mach-Zehnder geometry is realized with three material diffraction 
gratings. The total length of the interferometer including beam collimation is about 2.1 m and the molecules 

pass it at a mean velocity of 820 m/s. This yields an interrogation time of 2.6 ms. The maximally possible contrast is 
limited by two factors: First, the different weights of the interference paths, which correspond to the zeroth and the 
first diffraction order at the first grating; they are given by Pi/Pq = sine (0.37r) /l = 0.74. Second, the modulation of 
the interference pattern by the third grating mask; it contributes a factor of sine (0.37r) = 0.86 to the fringe amplitude, 
assuming a sinusoidal fringe pattern and a grating opening fraction of 30%. The detected contrast is thus limited to 
below 85%. We extract a measured fringe contrast of about 30% from the inset of Fig. 4 in [4(|, i.e. / ~ 0.35, which 
leads to fi = 7.2. 



Appendix E: Oscillating microresonators 



The authors of [30j propose to create a quantum superposition state of an oscillating micromirror by entangling it 
with a single cavity photon in one arm of a Michelson interferometer. If coherence is maintained in the mirror motion 
during one period of oscillation photon interference fringes should be observed at 100% contrast. 

The classicalization master equation (equation (1) in the main text) for the harmonic mirror motion can be inte- 
grated explicitly. We find that it reduces the fringe visibility after one oscillation period 2n/uj m by the factor 



R = exp 



2 71 



d£ 



0J m T 



9lD 



. 2 £ 2Hk . . 
2kxo sin — , sin £1—1 

2 x 



(El) 



with giD the reduced Fourier transform of the distribution g. The latter and the time parameter r are given by 
equations (6) and (7) in the main text. The micromirror is modeled as a homogeneous cube of mass density g = 



2300 kg/m 3 , b 



10 fxm edge length and a mass of M 

2 



gb 3 = 2.3 ng, which yields 



1 M 



7 



Sid (s,<?) = 1 
with the CTg-dependent geometry factor 



M 



7 



G q b 



2nm e cr s a l 



exp 



■ exp 



2h 2 



2m 2 a 2 



2a 2 



sine 



7T cr, 

2 



\V2h 



(E2) 
(E3) 

(E4) 



The authors presume a frequency uj rn /2ir = 500 Hz, a ground state oscillation amplitude of xq = 170 fm and a 
photon-mirror coupling strength of k — 1.63. We find /i = 19.0 for a measured 50% fidelity. 

For the hypothetical superposition experiment with an oscillating Al micromembrane, as listed in Table I of the 
main text, we use the parameters given in [47]. The membrane mass M = 48 pg and the mechanical frequency 
w m /27r = 10.56 MHz yield a tiny ground state amplitude of xq — \j2h/Muj m = 8 fm. To give a good upper estimate 
of the macroscopicity of such an experiment we thus approximate the flexural mode of the membrane by a axial 
center-of-mass vibration of a homogeneous disc of thickness b = 100 nm and radius R — 7.5 /im. We obtain the 
effective distribution <?id (s,<?) of (|E3|) and a time parameter 



^1 



M 



a q R 



qR 

h? 



exp 



<R 2 
h 2 



exp 



<R 2 



h 



h 2 



(E5) 



with Iq i the modified Bessel functions. 

We assume that the membrane is prepared in the superposition state l^) = (|0) + |1)) /\/2 of the zero- and the 
one-phonon eigenstate, and that the associated nondiagonal matrix element (l|pj|0) does not decay by more than 
50% after a time t = 27m/cj m which corresponds to n = 1000 oscillation cycles. Due to the large mass M and the 
small amplitude xq we may neglect the position spread a s in (|E3[) . and we may Taylor-expand its Fourier transform 
giu (x,p) to lowest order in x. With this we arrive at the explicit condition 



27TO 



(iHQ) 

(1|V)(V|0) U m r 7 6 2 



1 — exp 



2ft 2 



> 50%, 



(E6) 
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which leads to a macroscopicity of /i = 11.5. 



Appendix F: BEC interference 



The interference of two sodium BECs observed in Q is modeled using a second quantization phase-space picture 
l (52j . Following the same line we define a second quantization form of the characteristic function, 



dx e lpx 



(Fl) 



De Broglie interference of trapped BECs is observed as a fringe pattern in the time-evolved single-particle density 
n (x) = ijv (x) if} (x) for each individual run of the experiment. The fringe visibility is given by the corresponding 
Fourier component of n (x) , 



4 



dxi 



e 



2-rrixQ / A 



n(x ) . 



(F2) 



where A denotes the fringe spacing. The pattern observed in each run of the experiment [|[ can be assessed in the case 
of non-interacting bosons by replacing the annihilation operator ip (x) with the collective wave function ip (x) of the 
two trapped condensates. A free evolution of (|F1[) by the time t then yields the visibility \ (0, h/X) of the resulting 
interference pattern. The fringe spacing A = should be modified to account for interactions in the BEC |52j |. 
The second quantization form of the reduced TV-particle operators (4) reads as 



W («.«) = — fdxe- tqx/h iJ;Ux) 1 jj(x + — s) 
ra P / V m J 



A straightforward calculation reveals that (|F1[) then classicalizes at the rate 



£x(x,p) 



m 
m e 



1 

7V. 



ds dq g e 



-s,q 



i(qx-ps)/h 



(F3) 



(F4) 



of a single atom. We therefore estimate the macroscopicity from equation (9) of the main text. The authors of Q 
observed about / = 75% interference contrast in a sodium BEC after a time-of-flight of t = 40 ms; this yields /1 = 8.4. 
In the experiment 53] the phase sensitivity of the interferometer was increased, but at an interference contrast of 
only 15% after 200ms, which leads to /i = 8.3 (not discussed in the main text). 

Modern-day experiments with multi-component BECs make use of nonlinear interactions and number squeezing 
to increase the coherence time and the phase sensitivity employing internal atomic states [54T - I561 ] . The use of such 
techniques in interference experiments with spatially split BECs would only yield a macroscopicity fi that is comparable 
to single-atom interferometers. This is due to the fact that the single-particle nature of the classicalizing effect (IF4|) 
holds irrespectively of whether nonlinear interactions modify the coherent time evolution of the condensate wave 
function. Larger values of /1 could be achieved by increasing the fringe visibility and the time-of-flight in both 
single-atom and BEC experiments, possibly carried out in a microgravity environment. 



Appendix G: SQUID interference 

For the case of SQUID experiments we obtain the exclusion curve of the classicalization parameters in Fig. 1 
and the /Li-values in Fig. 2 and Tab. 1 by estimating the decay rate of a superposition state of macroscopically 
different supercurrents, i.e. different phases across the junctions in a Josephson loop. This was studied theoretically 
for spontaneous localization models in [32j . whose observable consequences are a special case of the classicalizing 
modification discussed here [23j . 

A state of finite current density \j) in a solid with electron density n e is described by a Fermi sphere, displaced by the 
momentum hkj, j — n e ehkj /m e . The undisplaced state |0) is taken to be the BCS ground state of the superconductor 
[13, [Hj]. It is characterized by the probability amplitudes Vf. {uk = \J\ — v'j,) of a Cooper pair (fe f, — k I) being 
occupied (unoccupied), 

h 2 - h 2 

a = -|l = kF J - (CD 

(fe 2 - k 2 F f + (2m e A k /h 2 ] 
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Here, hkp = m e VF denotes the Fermi momentum and A k the pairing energy. The latter is approximated in the usual 
way by the zero-temperature energy gap A = 1.76k B T c for electrons close to the Fermi level, \k 2 — k 2 F \ < 2m e uj D /h, 
and zero otherwise. The term ujp, denotes the Debye cutoff frequency of the material. We use the literature values 
k F = 1.18 At* A = 1.44 meV, Hlu d = 23.7 meV for Nb, and k F = 1.74 A -1 , A = 0.17 mcV, fuv D = 36.9 meV for Al, 
respectively (3lL l59l|. 

The second quantization form of the classicalization operators (equation (4) in the main text) for electrons reads as 

W e (s, fig) = Yl E ( fc ) a - ( k + l)- ( G2 ) 

The sum covers all discrete electron momentum states in a given volume V , each state occupying the elementary cell 
(27rfi,) 3 /V in momentum space. The classicalization kick distribution g e (s,q) must be discretized accordingly. We 
find that a superposition state p of two distinct supercurrents j 1 and j 2 decays at a rate 



dt(jMt)\3, 



<Jl|p(*)|j3 



= Tdiff + Tdcph (G3) 

t=o 



due to classicalization. This assumes that (ji\j 2 ) = 0) and that the net number SN of electrons occupying different 
states in each superposition branch [31| is large. The decay rate splits into two contributions. The first one is related 
to momentum diffusion, which requires that at least one elementary unit of momentum 2nh/V 1 ^ 3 is transferred. In 
the continuum limit J2 k ^A^ 71 ") 3 / we find 

r diff = (2^)3 r / [J d3sd3 9d 3 fc g e (s, hq) u k v k+q (u k v k+q + v k u k+q e i(2k+q)s ^ . (G4) 

? >7r/y 1 /3 

The expression is ultimately bounded by Tdiff < N/r e in the limit of arbitrarily strong momentum kicks, when all 
N = n e V conducting electrons can be transferred from one branch of the superposition to the other. It does not 
depend on the actual value of the supercurrents. 

The second contribution represents the dephasing that comes from the classicalization-induced position kicks, when 
no momentum redistribution of the electrons takes place, 



4V 2 H 3 

Tdcph = ^ 27r )6 r J J ^ S d ' iq 9e ( S ' h ^ ( X _ 6 * 



d 3 fc v 2 k e lk s 



(G5) 



Here, h5k = m e (j 1 — j 2 )/n e e denotes the difference in momentum displacement of the two current branches. It is 
orders of magnitude smaller than the Fermi momentum, and \8k ■ s\ -C 1 holds for any reasonable kick distribution 
g e . Hence the dephasing contribution scales quadratically with the net difference in occupation of the two displaced 
Fermi spheres, <5A^ = AN\Sk\/3kp- While this may potentially be significant for large SQUID geometries, the diffusion 
contribution dominates in all existing real-size experiments. 

Experimentally measured coherence times T 2 of such current superpositions provide an upper bound for the deca 
rate T. We estimate T 2 by the smallest observed frequency splitting in the experiments [l| (T 2 ~ 1 ns) and 
(T 2 w 10 ns); the authors of [6(| estimate T 2 « 15 ns. Classicalization parameters which lead to T > l/T 2 are then 
excluded by each experiment. This yields the SQUID curve in Fig. 1, as well as the /^-values plotted in Fig. 2; the latter 
are computed with the boundary condition a s < 1 A < H/a a , as discussed in the main text. The superconducting 
loop is spanned by L = 560 /zm of Nb in [l|, 20 /jm of Al in [6(j, and 180 /xm of Al in i2j. We assume the respective 
material cross sections as 5 /mi 2 , 36000 nm 2 , and 1 fim 2 . The experiment [6CJ yields a smaller macroscopicity, /x = 3.3, 
than [l| (jtx = 5.2) due to its smaller ring geometry. Only the greater value is included in Fig. 2. The large hypothetical 
SQUID in Tab. 1 of the main text is a 20 mm loop of 100 /1m 2 cross section with a coherence time of 1 ms. 

The actual values of the supercurrents do not influence the [i- values, since the dep hasing contribution is negligible 
in all cases. We use a current difference of I\ — I 2 = 3 fiA for Fig. 1, as given in [311 ] . 



Appendix H: Schrodinger's gedankenexperiment 

In our version of the famous gedankenexperiment, as listed in Tab. 1 of the main text, we consider the hypothetical 
superposition state of an ideal cat sitting at two places X\ and x 2 that are 10 cm apart. The center-of-mass coherence 
of the cat then decays like 

dt{x x \p\x 2 ) = z 



{Xl\p\X 2 ) T 



i Jd 3 sd 3 q g{s,q) (e lq (x2 - Xl)/h - (HI) 
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due to classicalization. We have neglected the weak position diffusion in the classicalization master equation (1) here. 
The mean coherence time of this state shall be 1 s. In order to evaluate the above decay rate using the r and the g 
of a compound, as defined by (6) and (7) in the main text, we model the cat as a homogeneous sphere of water with 
a mass of 4 kg. 
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